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WELLPOSEDNESS AND SINGULARITIES OF THE WATER WAVE 

EQUATIONS 

SIJUE WU 


Abstract. A class of water wave problems concerns the dynamics of the free interface 
separating an inviscid, incompressible and irrotational fluid, under the influence of grav¬ 
ity, from a zero-density region. In this note, we present some recent methods and ideas 
developed concerning the local and global wellposedness of these problems, the focus is 
on the structural aspect of the equations. 


1. Introduction 


( 1 . 1 ) 


A class of water wave problems concerns the motion of the interface separating an inviscid, 
incompressible, irrotational fluid, under the influence of gravity, from a region of zero density 
(i.e. air) in n-dimensional space. It is assumed that the fluid region is below the air region. 
Assume that the density of the fluid is 1, the gravitational field is — k, where k is the unit 
vector pointing in the upward vertical direction, and at time t > 0, the free interface is E(f), 
and the fluid occupies region SI (t). When surface tension is zero, the motion of the fluid is 
described by 

v* + (v • V)v = —k — VP on f l(t), t> 0, 
divv = 0, curlv = 0, on O(f), t> 0, 

P = 0, on E(f) 

(1, v) is tangent to the free surface (t, E(t)), 
where v is the fluid velocity, P is the fluid pressure. There is an important condition for 
these problems: 

dP 

(L2) 

pointwise on the interface, where n is the outward unit normal to the fluid interface E (t) 
[36j : it is well known that when surface tension is neglected and the Taylor sign condition 
EH fails, the water wave motion can be subject to the Taylor instability 

The study on water waves dates back centuries. Early mathematical works include Stokes 
(351 . Levi-Civita [28], and G.I. Taylor [36]. Nalimov [30], Yosihara [45] and Craig [14] proved 
local in time existence and uniqueness of solutions for the 2d water wave equation (11.11) for 
small initial data. In [39] 00], we showed that for dimensions n > 2, the strong Taylor sign 
condition 

dP 

- > c 0 > 0 (1.3) 

always holds for the infinite depth water wave problem EH), as long as the interface is non¬ 
self-intersecting and smooth; and the initial value problem of equation EH) is locally well- 
posed in Sobolev spaces H s , s > 4 for arbitrary given data. Since then, local wellposedness 
for water waves with additional effects such as the surface tension, bottom and non-zero 
vorticity, under the assumption EH0 were obtained (c.f. H EH HQ [23 EB [Ml 133121 ) ■ 

Alazard, Burq & Zuily ES proved local wellposedness of E3> in low regularity Sobolev 
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^When there is surface tension, or bottom, or vorticity, (II.3D does not hold, it needs to be assumed. 
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spaces where the interfaces are only in C 3 / 2 . Recently, we proved the almost global and 
global wellposedness of eh for small, smooth and localized initial data for dimensions n > 2 
mm-, Germain, Masmudi & Shatah obtained global existence for 3d water waves for a 
different class of small, smooth and localized data m- Our 2d almost global existence result 
has now been extended to global by Ionescu & Pusateri [22] and independently Alazard & 
Delort [3], see [18] EO] for an alternative proof. Finally, we mention our most recent work 
on two dimensional water waves with angled crests [ 43 J nn m, in which we showed that 
for water waves with angled crests, only the degenerate Taylor sign condition eh holds, 
with degeneracy at the singularities on the interface. We proved an a priori estimate [251 
and local existence [44j for the 2d water wave equation EH in this framework. 

The advances of water waves theory rely crucially on the understanding of the structure 
of the water wave equations. Indeed, eh is a nonlinear equation defined on moving do¬ 
mains, it is difficult to study it directly. A classical approach is to reduce from EH) to 
an equation on the interface, and study solutions of the interface equation. Then use the 
incompressibility and irrotationality of the velocity field to recover the velocity in the fluid 
domain by solving a boundary value problem for the Laplace equation. However the fluid 
interface equation is itself a fully nonlinear and nonlocal equation, its structure is not easy 
to understood. It is by achieving better understandings of the structure of this equation 
that has enabled us to apply analytical tools to deduce informations on the nature of the 
fluid motion. 

In this note we describe the approach in [391SD1 HD [42]. Our focus is on the structural 
aspect of the work. It is clear that the 2d case is structurally simpler than 3d. Our strategy 
has been to first understand the two dimensional case, taking advantage of complex analysis 
tools, in particular the Riemann mapping theorem, then use Clifford analysis to extend the 
2d results to 3d. We note that although Riemann mapping is not available in 3d, it is by 
using it that has enabled us to understand the 2d case well enough to develop an approach 
that extends to all dimensions. 

We consider solutions of the water wave equation EH in the setting where 

v(£,f) 0; as l£l °° 

and the interface £(f) tends to the horizontal plane at infinityH 

In section [2] we discuss the local wellposedness of EH, the focus is on deriving the 
quasilinear structure of equation eh, c.f. nil mi ng. In section [3j we consider the 
global in time behavior of solutions of EH in the regime of small waves, the focus is 
on understanding the nature of the nonlinearity of equation eh, c.f. mm~ We give 
some preparatory materials in the Appendices. In Appendix A, we give some basic analysis 
tools such as estimates for commutators and operators involved in the equations. These 
inequalities tell us how a certain term behaves in terms of estimates, whether it is of higher 
order, or lower order etc., and guide us in our derivations of the structure of the equation. 
In Appendix B we give some commutator identities that is used in our derivations. 

We use the following notations and conventions: [ A , B] := AB — BA is the commutator 
of operators A and B ; compositions are always in terms of the spatial variables and we write 
for / = /(-,*), g = g{-,t), := fog(-,t) := U g f{-,t ). 

2. Local wellposedness of the water wave equations 

From basic PDE theory we expect that in general, the Cauchy problem for a hyperbolic 
type PDE is locally solvable in Sobolev spaces, and we can solve it, for a short time, by 

2 The problem with velocity v(£,t) —>- (c, 0) as |£| —> oo can be reduced to one with v tends to 0 at 
infinity by the following observation: if (v, P) with E(f) : § = is a solution of 111.Il l, then 

V = v(£ + (c, 0)t, t) — c, P = P(f + (c, 0)t, t), with E(t) — (c, 0)t : £ = £(■, t) — ( c , 0)t 
is also a solution of 1 11.Il l. 
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energy estimates and an iterative argument; while the Cauchy problem for an elliptic type 
PDE is ill-posed. Hence in order to understand whether the water wave equation ed is 
uniquely solvable for a positive time period for arbitrary given Cauchy data, it is crucial to 
understand its quasi-linear structure. In this section, we derive the quasi-linear structure of 
the water wave equation ED, our focus is on the 2d case, since it is through a thorough 
understanding of this case that has enabled us to extend our work to 3d. We only give 
a brief description of how to extend the 2d derivations to 3d. We show that the strong 
Taylor inequality m always holds for C 1,7 interfaces, while for singular interfaces only 
the weak Taylor inequality ED holds; and this implies that the quasilinear structure of 
the water wave equation ED is of hyperbolic type in the regime of C 1,7 interfaces and 
respectively of degenerate hyperbolic type in the regime that includes singular free surfaces. 
The derivation given here is based on that in [39i 00]; due to the scope of this lecture note, 
we will only discuss the structural aspect of the work .39., 4Qj, and leave out the proof for 
the local wellposedness of ED- The interested reader may consult El HO] for the proof. 

2.1. The equation of the fluid interface in two space dimensions. In two space 
dimensions, we identify (x,y) with the complex number x + iy\ Rez, Im z are the real and 
imaginary parts of z; z = Rez — ilmz is the complex conjugate. 

Let the interface E(f) : z = z(a,t), a £ 1 be given by Lagrangian parameter a, so 
Zt(a,t) = v(z(a,f);£) is the velocity of the fluid particles on the interface, ztt{oi,i) = 
v t + (v • V)v(z(a, t); t) is the acceleration; notice that P = 0 on E(t) implies that VP is 
normal to E(t), therefore VP = —iaz a , where a = — p-| §17! we have from the first and 
third equation of ED that 

z tt +i = iaz a (2.1) 

The second equation of ED: div v = curl v = 0 implies that v is holomorphic in the fluid 
domain Lift)', hence z t is the boundary value of a holomorphic function in Ll(t). 

Let O C C be a domain with boundary E : z = z(a), a £ I, oriented clockwise. Let fj be 
the Hilbert transform associated to 9: 

£>/(o0 = — pv. [ / Z \^\a\ f(P)dP ( 2 - 2 ) 

m J z{a) - z(p) 

We have the following characterization of the trace of a holomorphic function on 9. 

Proposition 2.1. [2? a. Let g £ L p for some 1 < p < oo. Then g is the boundary value 
of a holomorphic function G on LI with G(z) —> 0 at infinity if and only if 

(/-$)<? = 0. (2.3) 

b. Let f £ L v for some 1 < p < oo. Then P nf '■= |(f + fj)/ is the boundary value of a 
holomorphic function © on LI, with ©(z) —> 0 as |z| —> oo. 

c. £1 = 0. 

Observe Proposition 12.11 gives that Sj 2 = I on L p . 

From Proposition 12.11 the second equation of ED is equivalent to Zt = Sjzt- Therefore 
the motion of the fluid interface E (t) : z = z(a,t) is given by0 

(z a + i = iaz a (24) 

[z t = Sjz t 

m is a fully nonlinear equation. To understand whether the equation is well posed, a usual 
strategy is to quasi-linearize the equation by differentiating. Notice that it can be hard to 
analyze the Hilbert transform S) in the second equation since it depends nonlinear nonlocally 
on the interface z = z(a,t), this motivates us to use the Riemann mapping (c.f. 1391)0 

^Equation (11.11) and equation (12.41) are equivalent, see |3911401 . 

4 In [30] [45], a quasilinear equation for the interface was derived for small and smooth waves in terms of 
the Lagrangian coordinates. In m, by using the Riemann mapping, a more concise quasilinear equation 




4 


SIJUE WU 


2.2. Riemann mappings and the quasi-linear structure of 2d water waves. Let 

<f>(-,i) : f l(t) —> P- be the Riemann mapping taking f l(t) to the lower half plane P_, 
satisfying Hindoo $ z (z,t) = 1. Let 

h(a;t) := $(z{a,t),t), 

so h : K. — > R is a homeomorphism. Let h" 1 be defined by 

h(h~ 1 (a', t), t) = a', a'el; 

and 


Z(a',t) := zo h Z t (a',t) := z t ° h 1 (a',t), Z tt (a', t) := z tt o h 


be the reparametrization of the position, velocity and acceleration of the interface in the 
Riemann mapping variable a!. Let 


Z i0l '(a',t) := d a >Z(a’,t), Z ti0t >(a',t) := d a <Z t {oi ,t), etc. 

We note that <f> _1 (a',i) = Z(a',t). Notice that v o <f> _1 : P_ — > C is holomorphic in 
the lower half plane P_ with vo = Z t (a',t). Precomposing (12.11) with h _1 and 

applying Proposition ^. II to vo$ _1 in P_, we have the free surface equation in the Riemann 
mapping variable: 

f ^tt + * = iAZ a i 

\z t = mz t 

where A o h = a h a and H is the Hilbert transform associated with the lower half plane P_: 


Uf(a') - -pv. [ —L—fWdp'. 

7TI J OC — p' 

0 As has been shown in [39], the quasi-linearization of (12.41) or (12.51) can be accomplished 
by just taking one derivative to t to equation m- 

Taking one derivative with respect to t to EH), we get 


zttt + iazt a = —ia t z a 


~(ztt-i), 

a 


( 2 . 6 ) 


the free surface equation is now (12.61) with the constraint Zt = Sjzt- Precomposing with h 1 
on both sides of (12.61) we have in the Riemann mapping variable the free surface equatior@ 


J Zttt + iAZt^i — ^ ° h 1 (Z t t — i) 

\z t = mzt 

We note that from the chain rule that for any function /, 


(2.7) 


U h -id t U h f(a',t) = (d t + Bd a ’)f(a',t) 

where B = h t o h -1 . Hence Z tt = (d t + Bd a >)Z t and Z ttt = (d t + Bd a ') 2 Z t . We will find 
B, A, and ^ o h~ 1 in terms of Z t and Z tt and show that m is a quasilinear equation for 


was derived for all non-self-intersecting waves. The approach described here is inspired by that in m , it 
was used when we extended the work in GS to 3d ESI- It is similar to that in [39] . with the difference that 
in m the derivation is in terms of the real component and here it is in both components. 

^The advantage of using the Riemann mapping is that in P_, the boundary value of a holomorphic 
function is characterized by g = Hp, see Proposition 1 2.1 j The kernel of the Hilbert transform HI is purely 
imaginary and is independent of the interface, it is easier to use HI to understand the relations among various 
quantities and hence the quasi-linear structure of the free surface equation. 

(12.61) and (12.71) are equivalent to d2.ll) and provided the initial data for d2.6f) and (|2.7[) stafisfy ED 
or H2.5D at t = 0. 
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U = Zt with the right hand side consisting of lower order terms 0 0 To this end, we need 
some basic estimates for commutators, we leave these and some other preparatory materials 
in the Appendices. The reader may want to consult the Appendices before continuing. 


2.2.1. The quantity A and the Taylor sign condition. The basic idea of deriving the formulas 
for A , B and h~ l is to use the holomorphicity or almost holomorphicity of our quantities, 
and the fact that Re(7—H)/ = / for real valued functions /. We will often write (7 — M)(fg) 
as a commutator: (7 — H)(/<?) = [/,H]g when g satisfies (7 — H )g = 0, since commutators 
are favorable in terms of estimates, see Appendix A. 

Let D a := j-d a , and D a i := -A—d a i. Notice that for any holomorphic function G in 

17(f) with boundary value g(a,t) = G(z(a, f), f), 

D a g(a,t ) = G z (z(a,t),t) 


and for any function /, 

(■ D a f ) o h" 1 = D a >(f o fi _1 ). 

We note that A is related to the important quantity — ^ by A o h = ah a and a = — ^, 

therefore 

= (2 8) 
In this subsection we derive a formula for the quantity A. This formula was first derived in 
139] to show that the strong Taylor sign condition m always holds for smooth non-self- 
intersecting interfaces. It has also played a key role in our recent work on 2d water waves 
with angled crests [25] [441. 

Taking complex conjugate of the first equation in (12.51) then multiply by Z a > yields 

Z^ a '{Ztt — i) = — iA\Z >a '\ 2 := —iA\. (2-9) 


The left hand side of (1^1) is almost holomorphic since Z >a / is the boundary value of the 
holomorphic function (ib” 1 ),/ and Zt.t is the time derivative of the holomorphic function 
Zt- We explore the almost holomorphicity of Ztt by expanding. Let F = v, we know F is 
holomorphic in 17(f), and Zt = F(z(a,t),t)> so 


ztt = F t (z(a, f), f) + F z (z(a,t),t)zt(u,t), z ta = F z (z{a,t),t)z a (cx,t) (2.10) 

therefore 

z tt = F t o z + — z t . (2-11) 

%Oi 

Precomposing with /i , subtracting —i, then multiplying by Z a /, we have 

Z.ot'i^tt ~ i) = Z (%' Ft O Z + ZtZt^a.' — iZ^a' = — iA\ 

Apply (7 — H) to both sides of the equation. Notice that F t oZ\s the boundary value of 
the holomorphic function F t o ih -1 , so (7 — M){Z^F t o Z) = 0, (7 — i)7 0 < = ljj therefore 

-i(7 - B.)A! = (7 - H ){Z t Z t , a >) - i 


7 Frorn Proposition 12.11 (/ — f))f measures how un-holomorphic a function / is. The reason (12~6] | is 
quasi-linear with the right hand side of lower order is that the two terms on the left hand side are ” almost 
holomorphic” in the sense that (/ — f))(zttt -\-iad a zt) is a commutator - since by (I — S))zt = 0, (I — + 

i&zta) = [dt + ivtdot-, $)\zt, while the conjugate of the right hand side points in the normal direction. We 
know a commutator is of lower order (c.f. Appendices A and B); and a holomorphic function with real part 
zero must be a constant; similarly an anti-holomorphic function with zero tangential part on the boundary 
must be a constant. If the left hand side of (| 2.6 1) were holomorphic, then it would have to be zero. Now 
the left hand side is almost holomorphic with (I — $))(zttt +iad a zt) a lower order term, then the right hand 
side must be a lower order term. We get this insight after our work in [39], not before. 

8 We sometimes abuse notation and say a function / is holomorphic if / is the boundary value of a 
holomorphic function in Q(£). 

^We know (<I > —1 ) z / —>■ 1 as z' —> oo; we assume a priori that (<I > —1 ) z rFt o$ —1 —» 0 as z' —> oo. As proved 
in m equation (l2.5l> is well posed in this regime. 
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Taking imaginary parts on both sides and using the fact (J — H )Zt, a ' = 00 to rewrite 
(/ - B.)(Z t Z t ,a') as [Z t ,W)Z t ,a> yields 


Ai = l-lm[Z t ,M\Z t ,c 


( 2 . 12 ) 


Compute 


- = -i Ref (W.0-W.<))<V(f.( dft 


i 


\Z,(a'A) - Z t (ii' 


a' — P' 

JL f \Z t (a',t)-Z t (p',t)\ 2 

27T ./ 


dp 


(2.13) 


dp' 


(a' — P') 2 

where in the last step we used integration by parts. We conclude 
Proposition 2.2. (c.f. [32]; Lemma 3.1) We have 1. 

1 r\Z t (a',t)-Z t (p',i)\ 2 


A\ — 1 — Im[Zt, — 1 + 


27T 


(o' - P') 2 


2 . 


dP A\ 


dn | Z a > j 


dP' > 1. (2.14) 


(2.15) 


in particular if the interface £(i) £ C 1,7 for some 7 > 0, then the strong Taylor sign 
condition <ra> holds. 


From (12.911 we have 

with Ai given by (12.141) . 


A = 


Ai 


I Z a - i\- 

A 1 


(2.16) 


2.2.2. The quantity B = ht o h 1 . The quantity B = ht o h 1 can be calculated similarly. 
Recall h(a,t ) = 4>(z(a, t), t), so 

h t = 4> t o z + (<F, o 2)z t , h a = (4> z o z)z a 

hence h t = 4> t 0 z + z Z t- Precomposing with h~ * l yields 

+ A (2.17) 

^ ,ot' 

Now (/ — H)$t o Z = (/ — H)(^-^ — 1) = 0 since 47 o Z and -A— — 1 are boundary values of 
the holomorphic functions 4> t o $ _1 and W — 1 respectively^ Apply (/ — H) to both 
sides of (12.171) then take the real parts; rewriting (/ — H )( Z t (-^-. _ 1)) as [ z t , H](^_ — 1), 
we get 

B = h t oh~ 1 = Re ( [Zt , H](—— - 1)) +2ReZ t . (2.18) 

^ ,a' 

Here we used the fact that Re(J — H )Z t = 2 R ,eZ t , since (/ + H )Z t = (I — M)Z t = 0. 


10 Because (I — H )Zt = 0. 

1 1 Again, we assume a priori that <t>t o <T<~ 1 (z'. t) — t 0 as 2 ' —> 00 . It has been proved in m that such 

solutions exist. 
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2.2.3. The quantity ^ ° h 1 . We analyze ^ ° h 1 similarly. Start with (12.61) or the first 
equation of m, and expand the left hand side using z t = F(z(a,t),t). Taking one more 
derivative to t to the first equation in (12.101) . we get 

zttt = ( F zz o z)z 2 + 2 (F tz o z)z t + (. F z o z)z tt + F tt o z. 

From (12.111) . the second equation of (12.101) . and the holomorphicity of F z , F t , we have 

F z ° z — D a Zf, F zz o z — D a z t , and F^ z o z — F) a (ztt [D a zt)zt) 


therefore 

Zttt = (D^zt)zt + 2 ztD a (ztt. — ( D a Zt)zt ) + ( D a zt)zu + F t t ° z. 

Precomposing with h~ x yields 

Zttt = (- D 2 a ,Zt)Z? + 2 Z t D a ,(Z tt - (D a iZ t )Z t ) + {D a ,Z t )Z tt + F tt o Z. (2.19) 

Now multiply the first equation of m by Z a i. then substitute in (12.191) . We have, by 
(12+1) . 

Z,a'{.{Da>Zt)Zt + 2 ZtD a i(Zt.t — {D a iZt)Zt) + {D a f Zt)Ztt+Fu o Z + iAZt : a'} 


= —iA\ — oh 1 . 


( 2 . 20 ) 


Apply (/ — H) to both sides of (12.201) . Using the fact that (I — H ){d a ’D a fZ t ) =0, (I — 
M)d a '(Z t t — (D a >Z t )Zt) = 0, (I — M)d a ‘Z t = 0 and (/ — H)F tt o Z = 0 because of the 
holomorphicity of d z >(F z o 'll -1 ), d z i(F t o $“ 1 ), d z r(F o $ _1 ) and F t t o 'll -1 , and the identity 
Ztt + * = iAZ ta i (12.51) . we rewrite each term on the left as commutator and get 

[Z 2 ,U]d a ,D a ,Z t + 2[Z u W]d a ,(Z tt - {D a ,Z t )Zt) + 2{Z tt ,U}d a ,Z t 

a, , (2.21) 

= (/-H)(-zAi— oh- 1 ). 

a 

For the sake of estimates we need to further rewrite the first term and the second part of 
the second term. Using integration by parts on each term, after cancelations we have 

1 [ (Z t (a',t)- Z t {p\t)) 2 


[Zf : U}d a ,D a ,Z t - 2[Z t ,W]d a ,((D a ,Z t )Z t ) = — 


(a' - /+) 2 


-DpZ t {p,t)dF 


•— ~[Zt, Z t ; D a fZt]. 


Therefore we can write (12.211) as 


at 


2[Z t ,M]Z tt:a , +2[Z tt ,U]d a ,Z t - [Z u Z t -D a ,Zt] = (/ - o h- 1 ). 

a 

Taking imaginary parts on both sides and dividing by —Ai yields 

Qt 0 ^-i _ ~ I m (2[Zt 7 ^\Ztt,a' + 2[Ztt,R]d a 'Z t — [Z t , Z t ;F) a /Z t ]) 
a A\ 


( 2 . 22 ) 


(2.23) 


2.2.4. The quasilinear equation. Sum up (12.71) . (12.141) . (12.161) . (12.181) and (12.231) . we have the 
equation for the free interface: 


where 


j (i % + Bd a ') 2 Z t + iAd a 'Zt — g 

\z t =MZ t 

Ztt = (dt + Bd a ,)Z t 

M = 1 - Im[Z t ,M\Z t , a , = 1+2 l -J W 

B = ht o h- 1 = Re ([Z t ,H](^i T - 1)) + 2ReZ t 


(2.24) 


(2.25) 


1 — • Z t t —i 

z a >~ M 


j = ( z tt - i) 


^-Im(2[Z t ,aiZ tta ,+2[Z tt ,Bld a ,Z t -[Z t ,Zf,D a ,Z t ]) 
A i 
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From the inequalities in Appendix A, we know in the regime where the interface is C * 1 ’ 7 , 
7 > 0, lj2.24M2.25ll is a quasilinear equation of the conjugate velocity Zt, with the right 
hand side of (12.24j) consisting of lower order terms. (12.24j) is of hyperbolic type since 


A = 



> cq > 0 


forC' 1 ’ 7 interfaces, and id a >Z t = |9 Q /|^0 Local well-posednessof (12.24l) - (12.25l) for ( Z t ,Z tt ) € 
C([0,T\,H s+1 / 2 x H* ), s > 4 has been proved by the energy method and an iterative argu¬ 
ment, we refer the reader to [35[ for details. 

We make the following remarks concerning some recent works. 


Remark 2.3. 1. The difference between the quasilinear equation (12.24M2.25D and the 
quasilinear equation (4.6)-(4.7) of [39] is that (12.24D - (12.25|) is in terms of both components 
of Z t and (4.6)-(4.7) of [31] is in terms of the real component of Z t . (4.6)-(4.7) of [31] has 
been written in a way so that it is easy to prove its equivalence with the interface equation 
(12.41) or equation (1.7)-(1.8) of [39]; see §6 of [39] . 

2. lj2.24M2.25ll is an equation for the conjugate velocity Z, t and conjugate acceleration 
Ztt, the interface doesn’t appear explicitly, so a solution of ll2.24M2.25jl can exist even when 
Z = Z(-,t ) becomes self- intersecting!*^ Checking through the derivation above we see that 
we arrived at (12.24|) - (12.25l) from (11.11) using only the following properties of the domain: 1. 
there is a conformal mapping taking the fluid region f2(f) to P_; 2. P = 0 on E(f). We note 
that ^ —> 2 1 / 2 is a conformal map that takes the region C \ {z = x + zO, x > 0} to the upper 
half plane; a domain with its boundary self-intersecting at the positive real axis can therefore 
be mapped conformally onto the lower half plane P_. Taking such a domain as the initial 
fluid domain, assuming P = 0 on E (t) even when E(f) self-intersects0 one can still solve 
equation (12.24M2.25p for a short time. Indeed this is the main idea in the work of [8]. Using 
this idea and the time reversibility of the water wave equation, by choosing an appropriate 
initial velocity field that pulls the initial domain apart, [8] proved the existence of ’’splash” 
and ’’splat” singularities starting from a smooth non-self-intersecting fluid region. 

3. The above derivation applies to fluid domains with arbitrary non-self-intersecting 
boundaries. We have from (j2.151) and (12.141) that the Taylor sign condition (|1.2D always holds, 
as long as the interface is non-self-intersecting. Assume that E(f) is non-self-intersecting 
with angled crests, assume the interior angle at a crest is v. We know the Riemann mapping 
4>~* (we move the singular point to the origin) behaves like 

~ (z'Y , with v = rn 


near the crest, so Z iCl > ss ( a') r 1 near the crest. From (USD and the fact A\ > 1 , we can 
conclude that the interior angle at the crest must be < n if the acceleration \Ztt\ ^ oo; 
we can also conclude that — ^ = 0 at the singularities where the interior angles are < n, 
therefore in the regime that includes singular free surfaces, the quasilinear equation (12.241) - 
(I2.25D is degenerate hyperbolic, c.f. [251135] . In [25] 33] we proved an a priori estimate and 
the local existence for 2d water waves in the regime including interfaces with angled crests, 
showing that the water wave equation DD admit such solutions. 

4. The quasi-linear equation for water waves in the periodic setting can be derived 
similarly, see [25] . 

5. The Riemann mapping variable is used in recent work |18 [ l20l . 


* 2 This is because id a i Zt = id a /MZt = \d a /\Zt. 

Z = Z(-,t ) is defined by z(-, t) = z(-, 0) + J 0 * z 3 (-, s) ds, where zt = Ztoh, Z = zofe -1 ; and ht = B(h, t). 

1 '\Ve note that when X)(t) self-intersects, the condition P = 0 on S(t) is unphysical. 
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2.3. The quasilinear equation in Lagrangian coordinates. In order to extend our 
work for 2d water waves to 3d, we need a derivation that does not rely on the Riemann 
mapping. Upon checking the derivation in d2.2l we see that all we have done is to apply 
Re(/ — H) to calculate the parameters in the equations. We certainly can do the same 
calculations with (/ — fj), Riemann mapping is not needed. 

We now use this idea to analyze cp| z a | and show that indeed (12.611 is quasilinear with the 
right hand side of lower order. 

First by EU and Proposition 12.21 we have 


o Z a = \Ztt + i\ 


and — i-. 


Ztt ~ i 


(2.26) 


\Za | | Ztt+l, 

We apply (/ — fj) to (12.61) . Using z t = Sjz t and Proposition IB. II in Appendix B, we have 
(I - Sj)(-ia t z a ) = {I - + iaz ta ) 

= [d t 2 +iad a ,$)\z t 

1 

iri 


= 2[ztt,fi\—— + 2[zt,fi}-—— 

Zry Zfy 




(2.27) 


z(a,t) - z(/3,t) 


Multiply both sides of (12.2711 by i and take the real parts. Since a and a t are real valued, 
we have 

(/+jr)(o t |2 a |) = 

1 _ (2.28) 


IZa 1 Zta r 1 Ztt.g J. I , Zt(a,t) Zt{/3,t) ^_ ,, 

R ^\z a \^ tU ^ z a Za 7T ij( z(a,t) - z(P,t) ) tftd ^ 


where 


8*f(a,t) = p.v. J Re{ 


-1 z n 




-}/G8,t)d/3 


ni \z a \ (z(a,t) - z(P,t)) 
is the adjoint of the double layer potential operator A in L 2 (E(f), dS). We know I + &* 
is invertible on L 2 (T,(t),dS) (cf. [T51 I23] h The second equation of (12.261) . (12.2811 and the 
estimates in Appendix A show that 0 t| 2 a | has the same regularity as that of Ztt and Zf 
We rewrite (EH) with the constraint zt = S)zt as 


Zttt 4“ i&Zta — \ztt+i\ ^ i II 
Zt = fiz t , 


(2.29) 


where a|z a | is given by (12.2611 and a t \z a \ is given by (12.28[) . Since z t is holomorphic, 
ijz—^d a zt = V n Z(. By Green’s identity, the Dirichlet-Neumann operator V n is a posi¬ 
tive operator. From Proposition 12.21 we know a\z a \ = —^ > cq > 0 in the regime 
of C 1 ’ 7 interfaces. Therefore (12.2911 is a quasilinear system of hyperbolic type, with the 
right hand side of the first equation in (12.291) consisting of terms of lower order deriva¬ 
tives of Zt- The local in time wellposedness of (12.29ll - (12.28l) - (12.26l) in Sobolev spaces (with 
(. Zt,Ztt ) G C([0,T],i7 s+1 / 2 x //(). s > 4 ) can then be proved by energy estimates and a 
fixed point iteration argument, [Jj 

This derivation has been extended to 3d using Clifford analysis, cf. [40], we will give a 
brief discussion on how to do this in ( 12.41 

Before ending this subsection, we mention that the quasilinear system (12.29l) - (12.28ll - (12.26ll 
is coordinate invariant. 

For fixed t, let k = k(a,t) : R -> R be a diffeomorphism with k a > 0. Let fc” 1 be such 
that k o A: -1 (a, t) = a. Define 

£ := z o/c -1 , b := kt o k _1 and Aok:=ak a . (2.30) 


lj A proof of the local wellposedness of the 3d counterpart of 112.2911 - 112.2811 - 112.26> is carried out in m- 
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Let 

D t := U^dtUk := d t + bd a (2.31) 

be the material derivative. By a simple application of the chain rule, we have 

u k 1( y d t + ™d a )U k = D 2 + iAd a , 

and equation (12.291) becomes 


\(D 2 +iAd a )D t (=(a t \z a \)ok 1 ^ + * 

(ac = hac 


(2.32) 


with 


and 


(I+IC*)((at\z a \) o fc _1 ) = 


Re&{2[Dfc H] + 2 [AC, H] 

I So: I S ol S« 


1 




A<(M) - D t C(P,t)s2 


C(<M) - COM) 




^ rr _i rfr ^ i / /(M)Cs(M) 
Hf(a,t) = U k SjU k f(a,t) = —p.v. / -=—-— -r—rdp, 

m J C{a,t) - C(/M) 

K* f{a, t) = p.v. [ Re{ — 

J t n 


IC^(/3M)| 


ICal (C(<M) - COM)) 


}f(P,t)d0 


(2.33) 


(2.34) 

(2.35) 


Notice the remarkable similarities between equations (12.29l) - (12.28l) and (12.32D - (12.33D . In 
particular, the structures of the terms in (I2.29I) - (I2.28D do not change under the change of 
variables. This makes it convenient for us to work in another coordinate system and to 
choose a different coordinate system when there is advantage to do so. In fact, this has been 
used in our study of the global in time behavior of water waves mm- 


2.4. The quasi-linear equation for 3d water waves. We derive the quasi-linear equa¬ 
tion for 3d water waves by carrying out the same procedure as for 2d. We first need to write 
down the 3d counterpart of the interface equation (12.41) . While equation (12.11) is readily 
available in 3d, to write down the second equation, we need a suitable counterpart in 3d of 
the equation for the trace on the interface of the velocity field v that satisfies div v = 0 and 
curlv = 0. This leads us to Clifford analysis. 

Let’s recall the basics of Clifford algebra, or in other words, the algebra of quaternions 
C(V 2 ) (c.f. ITT]). Let {I,ei,e 2 ,e 3 } be the basis of £(14), satisfying 

e 2 = — 1, fori = 1,2,3, e i e i = — e^e*, i^j, = eie 2 . (2.36) 

Let V = d x ei+d y e 2 + d z e 3 . By definition, a Clifford-valued function F : f2 C R 3 — > C(V 2 ) is 
Clifford analytic in domain f2 iff T>F = 0 in f2. Therefore F = JT =1 .fi e i is Clifford analytic 
in 12 if and only if divF = 0 and curlF = 0 in 12. Furthermore a function F is the trace of 
a Clifford analytic function in 12 if and only if F = fi^F, where 


tfsfl'O*, /?) = p.v. j J K(r)(a', 0) - 77 (a, $)) {rf a , x rf p ,)g{a', 0) da' d0 (2.37) 

is the 3d version of the Hilbert transform on E = <912 : 77 = r](a, (3), (a, 0) gR 2 , with normal 
Vc, x V /3 pointing out of 12, and 


m = 


1 

W3|r?| ’ 


2 77 

w 3 |77| 3 ’ 


K ( 77 ) = - 2 DT(t 7 ) 
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L 03 is the surface area of the unit sphere in R 3 0 As in the 2d case, if £ = £(ct,/3,t), 
(a, (3) £ R 2 is the free interface £(f) in Lagrangian coordinates (a, (3) at time t , with 
N = £ a x pointing out of the fluid domain, we can rewrite the 3D water wave system 
(11.11) (n = 3) as 


+ e 3 = aN 

6 = £>e (t)£t 


(2.38) 


where a= -^§£. 

Differentiating the first equation with respect to t yields 


ittt - aN t = a t N 
6 = 


(2.39) 


This is the 3d counterpart of the 2d quasilinear equation (12.61) with constraint Zt = F)Zt- It 
has been proved in [?D] that N t = — |JV|V n £t E and similar to the calculation (12.271) , (12.281) 
for 2d, an expression for a t \N\ in terms of (I-Sjs(t))(^ut+a\N\V a ^ t ) = [9 2 +a|A r |V n ,^ E ( t )]^ t 
has been derived. 

The Taylor sign condition ed also holds for C 2 interfaces of the water wave problem 
ed in 3d. This was proved in j4(Tj by an application of the Green’s identity. Here we give 
a heuristic argument via the maximum principle: 

Applying div to both sides of the Euler equation and using the assumption that curlv = 0 
yields 

AP = -|Vv| 2 <0 in n{t). 

Therefore from P = 0 on the interface £(£) and the maximum principle, we have — ^ > 0. 

(12.3911 is then a quasilinear equation of hyperbolic type with the right hand side consisting 
of lower order terms. The local in time wellposedness of (12.381) is proved in by applying 
energy estimates and an iterative argument to (12.391) . 


Remark 2.4. An analogous derivation of the quasilinear structure can now be performed 
directly on (O). Let D t = dt + v • V be the material derivative, we can rewrite the first 
equation in (ED as 

D t y + k = -VP. 

Apply D t to both sides to find 


P 2 v + [D t , V]P = —VD t P. (2.40) 

From the third equation: P = 0 on E(f), we have D t P = 0 on £(£), so VD t P points in the 
normal direction on the interface. We note that 

[D u d Xi ]P = -d x y ■ VP = -(VP • V) V j = a|A r |V n Vi 

on the fluid interface £(£)J3 so it corresponds to the term — aN t = alA^V^t in equation 
(12.39p . therefore (12.401) is the counterpart of our quasilinear equation (12.391) in the entire 
fluid region. In [lOl 129] . the authors used (12.401) to study a more general case where the 
vorticity curlv need not be zero FI 

We now turn to the question of long time behavior of solutions for the water wave equation 
ED for small initial data. 

16 (f2~37l) is similar to the 2d Hilbert transform EU in the sense that in 2d, the fundamental solution 
for Laplace equation is rfz) = In \z\, K(z) = —2(d x — idy ) (^: l n \ z \) — — the outward normal is iz a , 
and K(z(a) — z((3))(izp) = -K is the kernel of the Hilbert transform (| 2.2 1) . 

17 In 2d, we know iz a is a normal vector to S(t) pointing out of the fluid domain, and dt (iz a ) = —idaZt = 

12a | V nZt- 

In the second equality we used curlv = 0. 

■^In |101129| , the strong Taylor condition (11.3D is assumed to hold. 
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3. Global and almost global wellposedness of the water wave equations 

To understand the global in time behavior of the water wave motion, we need to under¬ 
stand the dispersion, the nature of the nonlinearity of the water waves and their interaction. 

In [HI in' we studied the water wave equation ED in two and three space dimensions for 
small data, we found a nonlinear transformation for the unknowns and a nonlinear change 
of the coordinates, so that the transformed quantities in the new coordinate system satisfy 
equations containing no quadratic nonlinear terms0 Using these canonical equations, we 
showed that for small, smooth and localized data of size e, the solution of the 2d water wave 
equation ED remain small and smooth for time 0 < t < e c / £ and for similar data, the 
solution of the 3d water wave equation ED remain small and smooth for all time. 

Let’s give a brief explanation of the dispersion of the water waves, the structural advan¬ 
tage of our canonical equations and how we found the transformations. 

Let u = z t (or u = £ t ). Linearizing the quasi-linear system (I2.29D - (I2.28D - (I2.26D (or (12.391) ') 
at the zero solution gives 

9 t 2 u + |D|u = F(u t , |D|u), (a,f)eP _1 xl (3.1) 

where |D| = y/—A, A is the Laplacian in R”^ 1 for n-dimensional water waves, F consists 
of the nonlinear terms. We know the dispersion relation of the linear water wave equation 

d 2 u+|-D|u = 0 (3.2) 

is 

w 2 = |fc| (3.3) 

for plane wave solution u = so waves of wave number k travel with phase velocity 

tc| k, where k = j|r, equation (13.11) is dispersive. For a large class *8 of smooth initial data, 
the solution of the linear equation (13.21) exists for all time and remains smooth, and its L°° 
norm decays with rate 1 /1 s . The question is for small data in for how long does the 
solution of the nonlinear equation ED remain smooth. We know nonlinear interactions can 
cause blow-up of the solutions at finite time. So to answer this question, we need to know 
for how long does the linear part of the equation (13.11) remain dominant. The weaker the 
nonlinear interaction, the longer the solution remains smooth. For small data, quadratic 
interactions are in general stronger than cubic and higher order interactions. 

To understand these assertions in qualitative terms, let’s consider the following model 
equation with a (p + l)th-order nonlinearity: 

<9 2 u + |D|u = ( d t u ) p+1 , (a, t) € R" _1 x R. (3.4) 

Suppose we can prove decay estimates for the solution: for i < s — 10, 

\d l d t u(t)\ L °o < (1 + E s (t) 1/2 , 

where d is some kind of derivatives, 

E s (t) = ^2 f \d J d t u(a,t)\ 2 + \d J \D\ 1/2 u(a,t)\ 2 da. 

Iil<s 

Then we can derive energy estimates for large enough s: 

jE a {t)<{i+t)- {J ^E s {ty/ 2+ \ 

therefore 

pT 

E S (Q)~ P ' 2 - E s {T)- p ' 2 < / {l + ty^^dt. (3.5) 

Jo 


2 ^When understood appropriately. 
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Heuristically, we expect to prove existence of solutions of (13.411 for as long as the energy 
E s (t) remains finite; by (13.511 . this can be achieved if 

{l + t)~ i ^dt <£ s (0)- p / 2 . 

Now if p = 1, i.e. if the nonlinear term in (13.4[1 is quadratic, then for both n = 2 and n = 3, 
the integral / 0 °° (1 + 1) 2 — 2 dt = 00 , and we would not be able to conclude solutions exist 
for all time for small initial data from this analysis. In fact for p = 1 and n = 2, the 
expected existence time is of order 0(e~ 2 ) for data of size e. If p > 2, i.e. if there is no 
quadratic nonlinearity, then we can expect to prove longer time existence for solutions of 
(13.41) for small initial data. In fact, we can expect for n = 2, p = 2 an existence time period 
of [0, e c / £ ]; and for n = 3, p = 2 an existence time period [0, 00 ) for data of size e, when e 
is sufficiently small. 

Now for the water wave equation m, the nonlinearity F(u t ,\D\u) contains quadratic 
terms, so appears too strong to conclude a global existence result. The question is whether 
there is another unknown t> that satisfies an equation of the type 

<9 t 2 0 + |D|0 = Fi(D t , |D|t>, «*, |£>|u) (3.6) 

with F\ containing no quadratic nonlinearities and ||o|| ~ ||u|| in various norms || • || involved 
in the analysis. The idea of finding such an unknown 0 is the so called method of normal 
forms, originally introduced by Poincare to solve ordinary differential equations. Certainly 
in most cases, one should not expect such a new unknown ti exist. For quadratic Klein- 
Gordon equation however, Simon and Shatah [321 f53] succeeded in finding a bilinear normal 
form transformation of the type 

0 = u + B(u, u) (3.7) 

with B( u,u) bilinear, canceling out the quadratic nonlinear terms in the Klein-Gordon 
equation, and satisfying the norm equivalence ||u + B(u, u)|| ~ ||u|| for small ||u||. 

For water wave equations (12.29l) - (12.28l) - (12.26l) (or (12.391) 1 however, a transformation of the 
type (13771) doesn’t quite work since it has a small divisor; working with the velocity potential 
and the Bernoulli equatioiQ a bilinear transform of the type (13.711 has a loss of derivatives 
(for detailed calculations and discussions, see Appendix C of [32].h What we did in pffl 132] 
was to further introduce a change of the coordinates. Indeed this makes sense since when 
one applies a method from ODE to PDE, it is reasonable to also take into consideration the 
spatial variables. 

We give in the next two subsections the transforms for the 2d and 3d water waves. The 
transforms are fully nonlinear. We found the transforms by first considering the 2d quasilin- 
ear equation (12.29H - (12.28D - (12.26D . starting with the ansatz (13.71) . taking into considerations 
of the coordinate invariance of (12.291) . looking for a coordinate invariant transformation! 2 ^ 
This entails much further efforts in understanding the bilinear transformation. What we 
finally arrived at is a fully nonlinear transform of the unknown function, coupled with a co¬ 
ordinate change. The process of finding the transforms is non-algorithmic. The transforms 
for the 3d water waves (12.391) is obtained by naturally extending the 2d version via Clifford 
analysis, c.f. 32] , 

3.1. The transformations for the 2D water waves. We give here the transformations 
we constructed in [44] for the 2D water waves. Let E (t) : z = z(a,t), agKbe the interface 
in Lagrangian coordinate a. 



21 Or equivalently the Zakharov-Craig-Sulem equation for the interface. 
22 The bilinear transform of the type (13.TI) is not coordinate invariant. 
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Proposition 3.1 (Proposition 2.3 of [41]). Let z = z(a, t) be a solution of the 2d water wave 
system (12.411 . Let II := (J— Sj)(z — z); let k = k(-, t) : M —> R be an arbitrary diffeomorphism. 
Let f := z o k- 1 , D t f '■= %t ° & _1 etc. be as in (12.30D - (12.31D . Then 

{D 2 t -iAdcfjiUok- 1 ) 

= *J -l«a,0-C(ft*)l a - t)d ° (3-8) 

2 r / D t ((a,t) -D t C(P,t) 

7T J V C(«,t) - C (P,t) 

Notice that the right hand side of equation (13.81) is cubicly small if the velocity D t f, the 
height function Im£ of the interface and their derivatives are small, but the left hand side of 
(13.81) still contains quadratic nonlinearities. Naturally we ask if there is a coordinate change 
k, so that b = k t o and A — 1 = (afc Q ) o fc -1 — 1 are quadratic. We need not look far, 
equation (12.181) suggests that we choose 

k(a, t) = 2 Rez(a, t) — h(a, t), a£l. (3-9) 


^ <9/3 Im C{P,t)d/3 


Proposition 3.2 (Proposition 2.4 of @4]). Let k be as given by (13.91) . b 
A = (a k a ) o k~ l . Let Tl be defined by (12.341) . We have 


(. I-H)b 
(I-H)(A- 1) 


So; 

i[ Dt <;,H]?^ + iiDtcn}^f± 

Set Set 


kt ok 1 and 


(3.10) 


Propositions 13.II and 13.21 show that the quantity 6 := II o k 1 = (/ — TL)(C ~ 0 with the 
coordinate change k given by (13.911 satisfies an equation of the type 

(df — id a )0 = Q 
(I + H)0 = 0 

with Q containing only nonlinear terms of cubic and higher order. We make the following 
remarks: 


Remark 3.3. 1. The transformation I — S) and the coordinate change k as given in (13.91) 
are fully nonlinear in terms of the unknown function z and its derivatives 0 

2. The bilinear part of the quantity II := (I — Sj)(z — z) has a bounded Fourier symbol. 
The coordinate change k takes care of the small divisor in the Fourier symbol of the bilinear 
normal form transformation. For a detailed explanation see Appendix C of [42] . 

3. For 9 satisfying (/ +H)9 = 0, (df — id a )9 = df9 + idaH0 = (df +1 D 1)0 + quadratic +.... 
The quadratic nonlinearity comes from H, which depends nonlinearly on the unknown if. 

4. Let (j) be the velocity potential, if>(a,t) := 4>{z{a,t),t) be the trace of <j) on the 
free interface. It has been shown in m that the quantities t/ fc -i(7 — Uk-idtJl also 
satisfy equations of the type (13.111) . with their equations given in Proposition 2.3 of [44] . 
daU^-i(I — SS)ip, U k -idtJi. l^aC/fc-iII and Imc^t/fc-iII are near identity transforms of the 
velocity D t f, 2D t (, the acceleration Dff and 2Im9 Q £, see Propositions 2.5 and 2.6 of |H] . 

5 . m-m is used in [38] to give a rigorous justification of the NLS from the 2d water 
wave equation (ED- 

6. The idea of changing coordinates is subsequently used in [13] to remove the quadratic 
nonlinear terms in the Burgers-Hilbert equation. 

We now give the proofs of Propositions 13.11 13.2pl 


“'’Tlia.l. is, 1 — U and k are not finite sums of multi-linear operators of z and its derivatives. 

2 ^The proofs are taken from m- 
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Proof. We first prove (13.8|) . Let z = z(-,t) be a solution of the water wave equation (12.411 . 
Apply (dt — iad a ) to II := (J — Sj)(z — z) and commute <9 2 — iad a with I — Sj to find 

(<9 t 2 - iad a ){(I - Sj)(z - z)} = (/ - /)){(<9 t 2 - iad a )(z - z)} - [<9 2 - iad a ,Sj](z - z). (3.12) 

Use (12.41) to find (9 2 — iad a )(z — z) = —2 z tt then use z t = Sjz t to write (/ — Sj)z t t as the 
commutator [d t . Sf\z t and applying Lemma \\>. 1 1 yields 

,Zta 


(I - i5){(5 2 - iad a )(z - z)} = (/ - Sj)(-2z u ) = -2 [z u Sj} : 
Applying Lemma IB. II to the second term gives 

i Zta Zta. 1 


(3.13) 


[df -iad ai f)\(z- z) =2 [z t ,$f\ : 


sZt{cK,t) Zt((3,t) 2 n / (n ,\ -in in 

*,j ( } 

(3.14) 


Subtract (13.141) from (13.131) . After cancelation this leaves, from (13.121) . 

(d?-iad a ){(I-f))(z-z)} = -2[z t ,Sj)^- + 1. j ( ^’^ _ ) 2 {z-z) 0 d/3. (3.15) 


z{a,t) - 

Because Zt and A — are holomorphic, from part 1 of Lemma 

z a 

[zt,fi-=-\Zta = = 0. 

%OL %OL 

In the first term of the right hand side of (13.151) insert [zt,f)f-\z ta to make it cubic. We 
have 

(<9 2 - ?a<9 a ){(/ - ff){z - z)} 


r 1 - 1 , 1 
— — 2[zt, •£)-h ]zta H-: 

7 Tl . 


,z t (a,t) - z t (P,t), i _ 


(3.16) 

□ 


Precomposing with k 1 and expanding the two terms on the right gives 
Proof. We prove Proposition 13.21 We have from (13.91) 

k — z = z — h. 

Recall h(a,t) = &(z(a,t),t) where <f>(-,f) : f2(f) -A P_ is the Riemann mapping satisfying 
lim-^oo $> z (z,t) = 1. We have 

h t = <L t o ^ + ($. o z)z t , h a = ($- o z)z a , 

therefore 

Zt - fa = $t o z + (& z o z - l)zt, z a - k a = (<E>- o z - l)z a . (3-17) 

Apply (/ — i^) to the first equality in (13.171) . Because <I> t , 4U are holomorphic in S2(f) 
with lim^^oo $> z (z,t) = 1, using Proposition 12.II and rewriting (/ — £)){( < iu o z — 1 )z t } as 
[z t ,Sj]($ z 02 - 1 ) yields 

-(/ - Sf)k t = (/ - S))(z t - k t ) 

= (-*"-■£>){($* °z~ = [zt,fi](®z 02 - 1 ) = [zt,fi] Za „ ka ■ 

Za 

Precomposing with fc -1 gives the first equality of (13.101) . 

Now multiply ia then apply (I — Sj) to the second equality in (13.171) . Using (12.41) and the 
fact that (/ — o z — 1) = 0, we also have 

(I-Sj)(iaz a - iak a ) = (/ - $))(iaz a ($ z o z - 1)) = [z u ,Sj](^ z oz - 1) = [z tt ,S )\———. 
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Use (12.411 and Lemma ft j. 1 1 to calculate 

(I - Sj)(iaz a ) = (I — S))(—z t t + i) = i- [d t ,Sj\z t = i- [z t ,Sj] — 

z a 

so 

-(I - Sj)(iak a ) = -(/ - Sj)(iaz a ) + [z tt ,Sj ]——— 

^ ( 3 - 1 9) 

= -* + [zt,$j]— + --• 

Z a Z a 

Precomposing with fc _1 yields the second equality of (13.101) . □ 

To extend the 2d coordinate change (13.91) to 3d, we need an expression that does not 
rely on the Riemann mapping. Observe that for the diffeomorphism k given by m, 
k — z = z — h and z — h is holomorphic with Im(z — h) = Imz, so we can replace z — h by 
\{I + Sj)(I + A)~ l {z — z), where A = Reij is the double layer potential operator, and 

k = z+ ^(I + S))(I + ^.)~ 1 (z — z) (3.20) 

modulo a real constant^ The expression (13.201) is directly extendable to 3d. 

3.2. The transformation for the 3D water waves. We extend the 2d transformations 
to 3d in the framework of the Clifford algebra C(V 2 ). Besides those in subsection 12.41 we 
need some additional notations. 

An element a £ C(V 2 ) can be represented uniquely by cr = ao + X)f =1 a i e ii with cq £ R for 
0 < i < 3. Define Re cr := ao and call it the real part of a. We call a a vector if ao = 0. If not 
specified, we always assume in an expression a = ao + Y^f=i th 0 ! cq £ M for 0 < i < 3. 
Define a := e 3 txe 3 , the conjugate of a. We identify a point or a vector £ = (aq, X 2 , y) £ R 3 
with its C(y 2 ) counterpart £ = aqei + + ye 3 . For vectors £, 77 £ C(V 2 ), we know 

& =-£-V + €*V, (3-21) 

where ^• 77 is the dot product, £ x r\ the cross product. For vectors £, r], £(£x ip is obtained 

by first finding the cross product £ x 77 , then regard it as a Clifford vector and calculating 
its multiplication with £ by the rule (12.361) . We write V = (d Xl ,d X2 ,d y ). We abbreviate 
notations such as 

f(a,/3) = JJ K(r)(a',P') - v(a,P)) iv' a ' x Vp) /(«', P') da'dp' 

:= JJ K( V ' - rj) (;q' a , x rj' pi ) f da'dp' := JJ K N' f da'd/3'. 

As in the 2d case, = / in L 2 , and ijsl = 0. 

We give the transformation for the 3D water wave equation (12.381) . Let the free interface 
£(t) be given by £ = £(a,/3,t) = 2 q(o,/?,f)ei + 2 : 2 ( 0 :, P,t ) e 2 + y(a,/3,t)e 3 in Lagrangian 
coordinates (o,/3) with N = x pointing out of the fluid domain Q(t). For fixed f, let 
k = k(-,t) = k\e\ + Aqe 2 : R 2 —>■ R 2 be a diffeomorphism with Jacobian J{k(t )) > 0. Let 
k~ x be such that k o fc _ 1 (o, /3, t ) = oei + / 3 e 2 - Dehne 

( — ^ok^ 1 , 6 -feoF 1 , A o kes := aJ{k)eo := ak a x kp (3.22) 

Let D t := U^dtUk be the material derivative, M := ( a x (p. By the chain rule, we know 

D t =d t + b-(d a ,d 0 ), Up 1 (aN x V)E4 = AM x V = A((pd a — ( a dp), (3.23) 


25 z — h and ^(/ + ft)(I + £) 1 (z — z) are holomorphic with the same imaginary part, so the difference 
between them is a constant in R. 
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and U k 1 S)-z(t)Uk ■= ?A(t), with 

n m f(a, P ,t) = JJ K(C(a', P', t) - ((a, p, t))(£, x CJg, )/(<*', A t) da'd^. (3.24) 

We have 

Proposition 3.4 (Proposition 1.3 of [42]). Let £ = £(a,/3,f) be a solution of the 3d water 
wave system (]2.38D . Let II = (/ — i5s(t))(C — P)> an d f or fixed t, k(-,t) : ffi 2 —> M 2 be a 
diffeomorphism. We have 

(Dt — AJ\f x V)(II o fc _1 ) 

= 2 JJ K(C - C) (AC - D'C) X {§,d a . - CAAAC 7 dc'd/3' 

- J [ K (C' - 0 (AC - ^CO X {{D' t C)f)>d a ' - (D(C%'^)(C' - CO da'd^ 

- JJ D t K(C - C) (AC - ^CO X (&3«' - Ca'A'XC' - CO da'd^ 

Observe that the second and third terms in the right hand side of (13.251) are cubicly 
small provided the velocity AC and the steepness of the height function d a ((— C), (C — C) 
are small, while the first term appears to be only quadratically small. Unlike the 2D case, 
multiplications of Clifford analytic functions are not necessarily analytic, so we cannot reduce 
the first term at the right hand side of equation (13.251) into a cubic form. However we note 
that the first term is almost analytic in the fluid domain f2(£), while the left hand side of 
(13.251) is almost analytic in the air region. The orthogonality of the projections \{I — 'HY,(t)) 
and |(/ + A(t)) allows us to reduce the first term into a cubic in energy estimates, see [4?j . 

Now the left hand side of (I3.25[) still contains quadratic terms. As in the 2D case, we 
resolve this difficulty by choosing an appropriate coordinate change k. Let 

k = k(a,P,t) = £(a,/3,t) - (I + io s(t) )y(a, P, t)e 3 + A s(t) i/(a, P, t)e 3 (3.26) 

Here A s(t) = ReJ5 E (t): 

•&£ (t)f(a,P,t) = - JJ K(£(a',P',t) - £{a,P,t)) ■ N'f(a , ,P',t)da' dp' (3.27) 

is the double layer potential operator. It is clear that the e 3 component of k as defined in 
(13.261) is zero. In addition, the real part of k is also zero. This is because 

J j m’ - £)x(& x ^,)y'e 3 da' dp’ 

= JJ (&& • K - • K)y’e 3 da' dp' 

= -2 JJ (&%T(£' - 0 - &0a'T(£' - 0)y'e 3 da' dp' 

= 2 JJ T(£' - rn'cVh' ~ Z'p'V^e* da' dp' 

= 2 JJ T(£' - OWei + N'e 2 ) da'dp' 

So 

^(t)ye 3 = %t)2/e 3 + 2 JJ T(£' - C)(7V( ei + N' 2 e 2 ) da' dp' (3.28) 

This shows that the mapping k defined in (13.261) has only the e± and e 2 components k = 
( ki,k 2 ) = kiei+k 2 e 2 . If E(t) is a graph with small steepness, i.e. if y a and yp are small, then 
the Jacobian of k = k(-,t ): J(k) = J(k(t)) = d a k\dpk 2 —d a k 2 dpk\ > 0 and k(-,t) : R 2 —> K 2 
defines a valid coordinate change (c.f. H2I). 
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The following proposition shows that if fc is as given in (13.261) . then b and A — 1 are 
quadratic. Let 

/C E (() P:=aei+/3e 2 , and ( ■—P + A. (3.29) 

Proposition 3.5 (Proposition 1.4 of [42]). Let k be as given in (13.261) . Let b = k t o fc -1 
and A o k = aJ(fc). We have 

b = \(H m - Wi^)AC -\[Du «E(t) - £s(t)](C - C) + ^s(t)(AC - AC) 

(A - l)e 3 = -(-%E(t) + T-L-e (t))L>?C + ^([-^’^swlAC _ [A,?A(t)]AC) 

+ x V,H s(t) ](C - C) - ^C/3 x (d Q /C s(t) (C - 0) 

+ ^^Ca x (dpJC-E(t)( C — 0) + ^A*A X C^A 

_ " (3.30) 

Are Hs(t) = ezH^t)^- 

Let % = II o fc -1 with fc be given by (13.261) . The left hand side of equation (13.251) is 

(df — e 2 (9 a + e\dp)x — dgXd a x + d a Xdpx + cubic and higher order terms 

The quadratic term dpXd a x — d a Xdpx is new in 3d. Observe that this is one of the null 
forms studied in [26]. It is also null for our equation and can be written as the factor 1/t 
times a quadratic expression involving some ’’invariant vector fields” for d 2 — e 2 <9 a + eidp, 
see [42]. Therefore this term is cubic in nature and equation (13.251) is of the type ”linear + 
cubic and higher order perturbations”. 

We refer the reader to [32] for the proof of Propositions 13.41 and 13.51 

We remark that the 3d transformations is recently used in Ezj to give a rigorous justi¬ 
fication of the modulation approximation for the 3d water wave equation m- Besides 
the 3d transforms, m uses the method of normal form to handle the quadratic term 
dgXd a x — d a Xdpx since this term is truly quadratic in the modulation regime. 

3.3. Global in time behavior of solutions for the 2d and 3d water waves. In 

[4TI 142] equations (|3.8I) and (|3.25l) together with the coordinate changes (13.91) and (13.261) 
are used to prove the almost global wellposedness of the 2d and global wellposedness of the 
3d water wave equations O) for small, smooth and localized initial data. The basic idea 
is what we illustrated with the model equation (13.41) : it is made rigorous by the method 
of invariant vector fields. This involves constructing invariant vector fields for the operator 
Of — e 2 d a + eidp (the invariant vector fields for d\ — id a for the 2d case is available due to the 
well studied Schrodinger operator idt — d%), proving generalized Sobolev inequalities that 
give L 2 —> L°° estimates with the decay rate l/t 1//2 for the 2d and 1/t for the 3d water waves, 
using equations (I3.8D - (13.91) and lj3.25M3.26ll to show that properly constructed energies that 
involve invariant vector fields remain bounded for the time period [0, e c / e ] for the 2d and for 
all time for the 3d water waves for data of size 0(e). The projection \{I —-£)) is used in various 
ways to project away ’’quadratic noises” in the course of deriving the energy estimates. We 
remark that it is more natural to treat D 2 — iAd a and D 2 — AN x V as the main operators 
for the 2d and 3d water wave equations than treating them as the perturbations of the linear 
operators d 2 — id a and d 2 — e 2 d a + e\dp. The almost global well-posedness for the 2d and 
global well-posedness for the 3d water wave equations follow from the local well-posedness 
results, the uniform boundedness of the energies and continuity arguments. For details of 
the proofs see Bum. We state the results. 

Let |D| = v 7 — A, FAR"- 1 ) = {/1 (/ + \D\)*f e L 2 (R n-1 )}, with ||/|| /r . (R »-i) = ||(J + 
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2d water waves. Let s > 12, max{[|] + 3,11} < l < s — 1. Assume 

z(a,0) = (a,y(a)), z t {a, 0) = u(a), z tt (a, 0) = ro(a) a e M, 

v (z, o) = g(z), zem (3 ' 31) 

and the data in (13.311) satisfies the 2d water wave system (12.41) . In particular g is a holo- 
morphic function in the initial fluid domain 12(0) and g(z(a,0)) = u(a). Let T = d a ,ad a . 
Assume that 

^ (l|r J ya||Hi/2(R) + l|r J u|| J/ 3 / 2 ( R ) + ||r J »|| ff i (R) ) < oo. 

bl<s-i 

Let 

e = X! (ll^yll^HR) + llr^ll-URR)) + X! ll("^) J glU 2 (f2(0))- 

\j\<l 3<l~ 2 

Theorem 3.6 (2d Theorem, c.f.ptTj). There exist eo and c > 0, such that for e < eg, the 
initial value problem (12.41) - (13.311) has a unique classical solution for the time period [0, e C//e ]. 
During this time, the interface is a graph , the solution is as regular as the initial data 
and remains small. Moreover the L°° norm of the steepness d a (z — z), the velocity Zt and 
acceleration Ztt decay at rate l/t 1 / 2 . 

3d water waves. Let s > 27, max{[|] + 1,17} < l < s — 10. Assume that initially 

£(<*,£, 0) = (a,P,y°(a,/3)), &(a, /?,0) = u°(a, /3), £ tt (a,P,0) = tu°(a,/3), (3.32) 

and the data in (13.321) satisfies the 3d water wave system (12.381) . Let T = d a , dp, ad a + 
Pdp, adp — Pd a . Assume that 

£ ||r J |i4| 1/2 y 0 || i 2( R 2) + ||r J 'd?/ 0 || jF/ i/ 2 ( R 2 ) + ||r J u 0 || jF/ 3/2 (R 2 ) + ||Pm 0 || ff i (R 2 ) < OO (3.33) 

d = d a ,dp 

Let 

e= £ ||r' 7 |D| 1 / 2 ?/ 0 || i 2 (R 2 ) + ||r J d2/°|| i 2 (R 2 ) + ||r J u°||// 1 /2 (R 2 ) + ||r J tD°||L2( R 2 ) . ( 3 . 34 ) 

lil<!+3 

B=a a ,Sp 

Theorem 3.7 (3d Theorem, c.f. [32]~). There exists eo > 0, such that for 0 < e < eg, the 
initial value problem (12.381) (13.321) has a unique classical solution globally in time. For each 
time 0 < t < oo, the interface is a graph, the solution has the same regularity as the initial 
data and remains small. Moreover the L°° norm of the steepness and the acceleration on 
the interface, the derivative of the velocity on the interface decay at rate \. 

Remark 3.8. 1. The existence time in Theorem 13.61 is extended to global in [2211512U] by 
further understanding the nature of the cubic nonlinearities in equation CD. 

2. In |41l a quick dispersive estimate was proved by the vector field method. An in- 
depth analysis on the dispersion of the linear water wave operator df + \D\ is performed in 
a threshold is found so that when the amount of small frequency waves in the initial 
data is below the threshold, the solution of the linear water wave equation m decays 
with rate 1/y/t , while above the threshold, there is a growth factor in the linear solution. 
Consequences in the nonlinear setting remain to be understood. 

Acknowledgement: The author thanks Jeffrey Rauch and Shuang Miao for carefully 
reading through the manuscript and for their remarks and suggestions. 
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Appendix A. Basic analysis preparations 


In this section we present some inequalities and identities on R that are used to guide the 
derivation of the quasilinear structure of the 2d water waves. Corresponding inequalities 
and identities are available in all dimensions R d , we refer the reader to [4(2 02] for those. 
Let H £ C 1 (R; R ; ), A, £ C^R), i = 1,.. .m, F £ C'°°(R i ). Define 


Ci{H, A, f)(x) =p.v. / F( 




x-y 


(i X - y ) m+1 


f(y) dy- 


(A.l) 


Proposition A.l. There exist constants Ci = ci(F,\\H'\\l°°), C2 = C2(F,\\H'\\l°°), such 
that 

1. For any f £ L 2 , A- £ L°°, 1 < i < m, 

\\Ci(H, A, /)|| L 2 <c 1 ||A' 1 |U-...||Ai ri || i -||/|| i 2 . (A.2) 

2. For any f £ L °°, A- £ L°° , 2 < i < to, A' x £ L 2 , 

l|C' 1 (H ) A,/)|| ia < C2 ||Ai|| La ||^|| L co...|| J C|| L o 0 ||/|| L «. (A.3) 


(IA.2I) is a result of Coifman, McIntosh and Meyer [T2] . (IA.3I) is a consequence of the Tb 
Theorem, a proof is given in US- 

Let U, A,. F satisfy the same assumptions as in (IA.1I) . Define 


C 2 (H,A,f)(x) 



H{x) - H(y) 
x-y 


HT=i(A,,(x) - A,{y)) 
(x - y) m 


d y f(y ) dy. 


(A.4) 


We have the following inequalities. 


Proposition A. 2 . There exist constants C3 = c 3 (F,\\F['\\ L c X ,), C4 = C4(P, such 

that 

1. For any f £ L 2 , A' £ L°°, 1 < i < m, 

||C 2 (//,A,/)|| i2 < caKIUoo... IK.IU-II/IU- ( A -5) 

2. For any f £ L°°, A[ £ L°° , 2 < i < m, A[ £ L 2 , 

||C' 2 (ff,A,/)|| L2 <c 4 ||A' 1 || i2 ||A'||^...||A;j|^||/|| L ». (A.6) 


Using integration by parts, the operator C^IL, A, /) can be easily converted into a sum 
of operators of the form Ci(H,A,f). (IA.5I) and (IA.6I) follow from (IA.2I) and (IA.3I) . 

The following identities are useful to compute the derivatives of the integral operators. 
Let 

K f(x, t ) = p.v. J K(,x , y, t)f{y , t) dy 

where either K or {x — y)K(x , y\ t ) is continuous and bounded, and K is smooth away from 
the diagonal A = {(a;, y)\x = y}. We have for / £ C' 1 (I? 1+1 ) vanishing as |x| —> 00 , 


[d t ,K]f(x,t) = J d t K(x,y;t)f{y,t)dy 
[d x ,K\f(x,t) = J(d x + d y )K(x, y\ t)f(y, t) dy 


(A.7) 


The first identity in (1A.7I) is straightforward, the second is obtained by integration by 
parts. 

Using (IA.7I) . 9®,[/, W[d a >g equals to the sum of [9®,/,H]9 a /gl3 [/,IHt] 1 r? and some 
intermediate terms. By Propositions IA.1I and IA.2I and the Sobolev embedding, [/, M\d a i<j 
has the same regularity as / and g. 


2 ®Expand \d^, f,W]d a /g = d®,/H(d a /g) — H (d^, fd a r g) and estimate term by term. 
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When the Riemann mapping is not used, typically we work in the regime where the 
interface is chord-arc, that is there are constants /ii > 0 and /12 > 0, such that 

/ii|a - P\< | z(a, t) - z(p, t)\ < fj, 2 \a - P\, Va, /? € R. 

By (1A.7[) . derivatives of [/, #]%£ are of types (IA.1I) and (1A.4[) . where in (IA.1I) or (IA.4I) 
F £ Cq°(C) is chosen to be F{z) = for /xi < |^| < n 2 and if (a) = z(a,t). An 
application of Proposition IA.1I and IA.2I and Sobolev embedding shows that [/, io] has 
the same regularity as f, g and z. 


Appendix B. Commutator identities 


We need the following identities in our derivation of the quasilinear structure (12.2711 and 
Propositions 13.11 and 13.21 Let z = z(-,t) define a non-selfintersecting curve, and be the 
Hilbert transform as defined in (12.211 . 


Lemma B.l (Lemma 2.1 of [51]). Assume that zt , z a — 1 £ C' 1 ([0, T], H 1 ( R.)), / G C ll (M x 
(0, T)) satisfies f a (a,t) —> 0, as |a| —> 00 . We have 


[dt,S)]f = [z t ,Sj] 


fa 


m2 f i r r r i/a 1 or fta 1 f /2j(c¥, t) Zt((3,t)\ 2„ 

[ 3 , ,W = [*.,«]- + «!---/ ( , (a , 0 _ zW0 ) * 




(B.l) 


[d 2 t - iad a ,$)\f = 2 [z t ,Sj\^ - — f d p 

Z a 7TI J ^ +\ ' 


Zto 


(/ - Sj)(-ia t z a ) = 2 [z tt ,Sj] -h 2[zt,fj] 


z(a,t) - z(fi,t) 

1 Z-tt 


1 

iri 


Zt{ot,t) - z t {p,t)^ 2 


z(a,t) - z(/3,t) 


) Zt/3 dp 


The proof of Lemma IB.II is straightforward by integration by parts. We omit the proof. 
Let fi C C be a domain with boundary T. \ z = z(a),a £ I oriented clockwise. Let S) be 
defined by (12.21) . 


Lemma B.2 (Lemma 2.2 of [41]). 1. If f = $jf , g = S)g, then [f,fj\g = 0. 
2. For any fig £ L 2 , we have [f,Sj]fig = ~[Sjf,Sj]g. 


The first statement follows from the fact that the product of holomorphic functions is 
holomorphic. Observe that the first statement also holds for / = —S)f and g = — the 
second statement follows from the first by applying the first identity to (/d=f))/ and (I±S))g. 

We do not give the commutator identities for 3d, but refer the reader to [40l SI] for 
details. 
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